Continuum Limit Research Articles

On the continuum limit of epidemiological models on graphs: convergence and approximation results

We focus on an epidemiological model (the archetypical SIR system) defined on graphs and study the asymptotic behavior of the solutions as the number of vertices in the graph diverges. By relying on the theory of graphons we provide a characterization of the limit and establish convergence results. We also provide approximation results for both deterministic and random discretizations.

Intrinsic rotation modulation by diffusive neutral particles in tokamaks

The modulated transport model, a model kinetic ion transport equation for the pedestal and scrape-off layer (SOL), is generalized to self-consistently include effects of a single neutral particle species. The neutrals contribute additional transport terms, modifying the v∥ -dependent orbit-averaged ion diffusivities of the original work and the resulting predicted intrinsic rotation of the ions. After making simplifying assumptions of the neutral transport, in particular taking the continuous transport limit via a short charge-exchange step expansion, we derive relatively simple analytic expressions that capture the diffusive neutral physics. Within the scope of the model’s validity, the neutral-driven intrinsic rotation can compete with the turbulence-driven intrinsic rotation. However, for physically motivated parameters, the neutral-driven intrinsic rotation appears negligible, either on a term-by-term basis or due to a strong cancellation between the neutral-driven momentum diffusion and pinch terms. It appears that a treatment containing finite charge-exchange steps is necessary to capture neutral transport of strong flow momentum into the confined region from the SOL.

Lattice regularization of reduced Kähler-Dirac fermions and connections to chiral fermions

We show how a path integral for reduced Kähler-Dirac fermions suffers from a phase ambiguity associated with the fermion measure that is an analog of the measure problem seen for chiral fermions. However, unlike the case of chiral symmetry, a doubler free lattice action exists which is invariant under the corresponding onsite symmetry. This allows for a clear diagnosis and solution to the problem using mirror fermions resulting in a unique gauge invariant measure. By introducing an appropriate set of Yukawa interactions which are consistent with ’t Hooft anomaly cancellation we conjecture the mirrors can be decoupled from low energy physics. Moreover, the minimal such Kähler-Dirac mirror model yields a light sector which corresponds, in the flat space continuum limit, to the Pati-Salam GUT model.

Quantum Geometrodynamics Revived I. Classical Constraint Algebra

Abstract In this series of papers, we present a set of methods to revive quantum geometrodynamics which encountered numerous mathematical and conceptual challenges in its original form promoted by Wheeler and De Witt. In this paper, we introduce the regularization scheme on which we base the subsequent quantization and continuum limit of the theory. Specifically, we employ the set of piecewise constant fields as the phase space of classical geometrodynamics, resulting in a theory with finitely many degrees of freedom of the spatial metric field. As this representation effectively corresponds to a lattice theory, we can utilize well–known techniques to depict the constraints and their algebra on the lattice. We are able to compute the lattice corrections to the constraint algebra. This model can now be quantized using the usual methods of finite–dimensional quantum mechanics, as we demonstrate in the following paper. The application of the continuum limit is the subject of a future publication.

Cosmology in Lorentzian Regge calculus: causality violations, massless scalar field and discrete dynamics

We develop a model of spatially flat, homogeneous and isotropic cosmology in Lorentzian Regge calculus, employing four-dimensional Lorentzian frusta as building blocks. By examining the causal structure of the discrete spacetimes obtained by gluing such four-frusta in spatial and temporal direction, we find causality violations if the sub-cells connecting spatial slices are spacelike. A Wick rotation to the Euclidean theory can be defined globally by a complexification of the variables and an analytic continuation of the action. Introducing a discrete free massless scalar field, we study its equations of motion and show that it evolves monotonically. Furthermore, in a continuum limit, we obtain the equations of a homogeneous scalar field on a spatially flat Friedmann background. Vacuum solutions to the causally regular Regge equations are static and flat and show a restoration of time reparametrisation invariance. In the presence of a scalar field, the height of a frustum is a dynamical variable that has a solution if causality violations are absent and if an inequality relating geometric and matter boundary data is satisfied. Edge lengths of cubes evolve monotonically, yielding a contracting or an expanding branch of the Universe. In a small deficit angle expansion, the system can be deparametrised via the scalar field and a continuum limit of the discrete theory can be defined which we show to yield the relational Friedmann equation. These properties are obstructed if higher orders of the deficit angle are taken into account. Our results suggest that the inclusion of timelike sub-cells is necessary for a causally regular classical evolution in this symmetry restricted setting. Ultimately, this works serves as a basis for forthcoming investigations on the cosmological path integral within the framework of effective spin foams.

Continuum Limit of the Green Function in Scaled Affine φ44 Quantum Euclidean Covariant Relativistic Field Theory

Through path integral Monte Carlo computer experiments, we prove that the affine quantization of the φ44-scaled Euclidean covariant relativistic scalar field theory is a valid quantum field theory with a well-defined continuum limit of the one- and two-point functions. Affine quantization leads to a completely satisfactory quantization of field theories in situations involving scaled behavior, leading to an unexpected term, ℏ2/φ2, which arises only in the quantum aspects.

A wealth distribution model with a non-Maxwellian collision kernel

Abstract A non-Maxwellian collision kernel is employed to study the evolution of wealth distribution in a multiagent society. The collision kernel divides agents into two different groups under certain conditions. Applying the kinetic theory of rarefied gases, we construct a two-group kinetic model for the evolution of wealth distribution. Under the continuous trading limit, the Fokker-Planck equation is derived and its steady-state solution is obtained. For the non-Maxwellian collision kernel, we find a suitable redistribution operator to match the taxation. Our results illustrate that taxation and redistribution have the property to change the Pareto index.

Connection between the symmetric discrete AKP system and bilinear ABS lattice equations

In this paper, we show that all the bilinear Adler–Bobenko–Suris (ABS) equations (except Q2 and Q4) can be obtained from symmetric discrete AKP system by taking proper reductions and continuum limits. Among the bilinear ABS equations, a simpler bilinear form of the ABS H2 equation is given. In addition, an 8-point 3-dimensional lattice equation and an 8-point 4-dimensional lattice equation are obtained as by-products. Both of them can be considered as extensions of the symmetric discrete AKP equation.

Large-N Principal Chiral Model in Arbitrary External Fields.

We report the explicit solution for the vacuum state of the two-dimensional SU(N) principal chiral model at large N for an arbitrary set of chemical potentials and any interaction strength, a unique result of such kind for an asymptotically free quantum field theory. The solution matches one-loop perturbative calculation at weak coupling, and in the opposite strong-coupling regime exhibits an emergent spacial dimension from the continuum limit of the SU(N) Dynkin diagram.

New solutions for the generalized q-deformed wave equation with q-translation symmetry

Abstract In this work, we explore the generalized discrete wave equation, which utilizes a specific irregular space interval. The introduction of this irregular space interval is motivated by its connection to the q-addition, a mathematical operation that arises in the nonextensive entropy theory. By taking the continuous limit, we obtain the wave equation with q-deformation, which captures the effects of the q-addition. To solve the generalized q-deformed wave equation, we investigate three different methods: the separation method, the reduced differential transform method, and the finite difference method. These methods offer distinct approaches for finding solutions to the equation. By comparing the results obtained from each method, we can evaluate their effectiveness and identify their respective strengths and limitations in solving the generalized q-deformed wave equation. The solutions obtained from this newly defined equation have potential applications in modeling physical systems with violated symmetries. The inclusion of the q-deformation allows for a more comprehensive description of such systems, which may exhibit nonextensive behavior or possess irregularities in their spatial intervals. By incorporating these features into the wave equation, we can improve our understanding and modeling capabilities of complex physical phenomena.

Beyond braid statistics: Constructing a lattice model for anyons with exchange statistics intrinsic to one dimension

Anyons obeying fractional exchange statistics arise naturally in two dimensions: Hard-core two-body constraints make the configuration space of particles not simply-connected. The braid group describes how topologically-inequivalent exchange paths can be associated to non-trivial geometric phases for abelian anyons. Braid-anyon exchange statistics can also be found in one dimension (1D), but this requires broken Galilean invariance to distinguish different ways for two anyons to exchange. However, recently it was shown that an alternative form of exchange statistics can occur in 1D because hard-core three-body constraints also make the configuration space not simply-connected. Instead of the braid group, the topology of exchange paths and their associated non-trivial geometric phases are described by the traid group. In this article we propose a first concrete model realizing this alternative form of anyonic exchange statistics. Starting from a bosonic lattice model that implements the desired geometric phases with number-dependent Peierls phases, we then define anyonic operators so that the kinetic energy term in the Hamiltonian becomes local and quadratic with respect to them. The ground-state of this traid-anyon-Hubbard model exhibits several indications of exchange statistics intermediate between bosons and fermions, as well as signs of emergent approximate Haldane exclusion statistics. The continuum limit results in a Galilean invariant Hamiltonian with eigenstates that correspond to previously constructed continuum wave functions for traid anyons. This provides not only an a-posteriori justification of our lattice model, but also shows that our construction serves as an intuitive approach to traid anyons, i.e. anyons intrinsic to 1D.

Continuum limit of the mobility edge and taste-degeneracy effects in high-temperature lattice QCD with staggered quarks

Continuum limit of the mobility edge and taste-degeneracy effects in high-temperature lattice QCD with staggered quarks

Immiscible Two-Phase Flow in Porous Media: Effective Rheology in the Continuum Limit

AbstractWe consider steady-state immiscible and incompressible two-phase flow in porous media. It is becoming increasingly clear that there is a flow regime where the volumetric flow rate depends on the pressure gradient as a power law with an exponent larger than one. This occurs when the capillary forces and viscous forces compete. At higher flow rates, where the viscous forces dominate, the volumetric flow rate depends linearly on the pressure gradient. This means that there is a crossover pressure gradient that separates these two flow regimes. At small enough pressure gradient, the capillary forces dominate. If one or both of the immiscible fluids percolate, the volumetric flow rate will then depend linearly on the pressure gradient as the interfaces will not move. If none of the fluids percolate, there will be a minimum pressure gradient threshold to mobilize the interfaces and thereby get the fluids moving. We now imagine a core sample of a given size. The question we pose is what happens to the crossover pressure gradient that separates the power-law regime from the high-flow rate linear regime and the threshold pressure gradient that blocks the flow at low pressure gradients when the size of the core sample is increased. Based on analytical calculations using the capillary bundle model and on numerical simulations using a dynamical pore-network model, we find that the crossover pressure gradient and the threshold pressure gradient decrease with two distinct power laws in the size. This means that the power-law regime disappears in the continuum limit where the pores are infinitely small compared to the sample size.

Cellular automaton ontology, bits, qubits and the Dirac equation

Cornerstones of the Cellular Automaton Interpretation of Quantum Mechanics are its ontological states that evolve by permutations, in this way never creating would-be quantum mechanical superposition states. We review and illustrate this with a classical Ising spin chain. It is shown that it can be related to the Weyl equation in the continuum limit. Yet, the model of discrete spins or bits unavoidably becomes a model of qubits by generating superpositions, if only slightly deformed. We study modifications of its signal velocity which, however, do not relate to mass terms. To incorporate the latter, we consider the Dirac equation in 1+1 dimensions and sketch an underlying discrete deterministic “necklace of necklaces” automaton that qualifies as ontological.

Majorana chain and Ising model - (non-invertible) translations, anomalies, and emanant symmetries

We study the symmetries of closed Majorana chains in 1+1d, including the translation, fermion parity, spatial parity, and time-reversal symmetries. The algebra of the symmetry operators is realized projectively on the Hilbert space, signaling anomalies on the lattice, and constraining the long-distance behavior. In the special case of the free Hamiltonian (and small deformations thereof), the continuum limit is the 1+1d free Majorana CFT. Its continuum chiral fermion parity (-1)^{F_L}(−1)FL emanates from the lattice translation symmetry. We find a lattice precursor of its mod 8 ’t Hooft anomaly. Using a Jordan-Wigner transformation, we sum over the spin structures of the lattice model (a procedure known as the GSO projection), while carefully tracking the global symmetries. In the resulting bosonic model of Ising spins, the Majorana translation operator leads to a non-invertible lattice translation symmetry at the critical point. The non-invertible Kramers-Wannier duality operator of the continuum Ising CFT emanates from this non-invertible lattice translation of the transverse-field Ising model.

Opinion Models, Election Data, and Political Theory.

A unifying setup for opinion models originating in statistical physics and stochastic opinion dynamics are developed and used to analyze election data. The results are interpreted in the light of political theory. We investigate the connection between Potts (Curie-Weiss) models and stochastic opinion models in the view of the Boltzmann distribution and stochastic Glauber dynamics. We particularly find that the q-voter model can be considered as a natural extension of the Zealot model, which is adapted by Lagrangian parameters. We also discuss weak and strong effects (also called extensive and nonextensive) continuum limits for the models. The results are used to compare the Curie-Weiss model, two q-voter models (weak and strong effects), and a reinforcement model (weak effects) in explaining electoral outcomes in four western democracies (United States, Great Britain, France, and Germany). We find that particularly the weak effects models are able to fit the data (Kolmogorov-Smirnov test) where the weak effects reinforcement model performs best (AIC). Additionally, we show how the institutional structure shapes the process of opinion formation. By focusing on the dynamics of opinion formation preceding the act of voting, the models discussed in this paper give insights both into the empirical explanation of elections as such, as well as important aspects of the theory of democracy. Therefore, this paper shows the usefulness of an interdisciplinary approach in studying real world political outcomes by using mathematical models.

Noise-induced collective actuation in active solids.

Collective actuation describes the spontaneous synchronized oscillations taking place in active solids when the elasto-active feedback, which generically couples the reorientation of the active forces and the elastic stress, is large enough. In the absence of noise, collective actuation takes the form of a strong condensation of the dynamics on a specific pair of modes and their generalized harmonics. Here we report experiments conducted with centimetric active elastic structures, where collective oscillation takes place along the single lowest energy mode of the system, gapped from the other modes because of the system's geometry. Combining the numerical and theoretical analysis of an agent-based model, we demonstrate that this form of collective actuation is noise-induced. The effect of the noise is first analyzed in a single-particle toy model that reveals the interplay between the noise and the specific structure of the phase space. We then show that in the continuous limit, any finite amount of noise turns this new form of transition to collective actuation into a bona fide supercritical Hopf bifurcation.

The chiral separation effect from lattice QCD at the physical point

In this paper we study the Chiral Separation Effect by means of first-principles lattice QCD simulations. For the first time in the literature, we determine the continuum limit of the associated conductivity using 2+1 flavors of dynamical staggered quarks at physical masses. The results reveal a suppression of the conductivity in the confined phase and a gradual enhancement toward the perturbative value for high temperatures. In addition to our dynamical setup, we also investigate the impact of the quenched approximation on the conductivity, using both staggered and Wilson quarks. Finally, we highlight the relevance of employing conserved vector and anomalous axial currents in the lattice simulations.

Hamiltonian and recursion operators for a discrete analogue of the Kaup-Kupershmidt equation

In this paper we study the algebraic properties of a new integrable differential-difference equation. This equation can be seen as a deformation of the modified Narita-Itoh-Bogoyavlensky equation and has the Kaup-Kupershmidt equation in its continuous limit. Using its Lax representation we explicitly construct a recursion operator for this equation and prove that it is a Nijenhuis operator. Moreover, we present the bi-Hamiltonian structures for this new equation.

Remarks on discrete Dirac operators and their continuum limits

Remarks on discrete Dirac operators and their continuum limits